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A  COMPARATIVE  STUDY  OF  TWO  AZIMUTH-BASED  NON-STANDARD  LOCATION  METHODS 

R.  Jih 

Department  of  State  /Arms  Control,  Verification,  and  Compliance  Bureau,  Washington,  U.S.A. 

In  this  study,  two  non-standard,  independently  developed,  azimuth-based  seismic  triangulation  procedures  are 
implemented  and  evaluated.  The  goal  is  to  develop  ad  hoc  algorithms  which  can  effectively  tackle  the  seismic  location 
problem  under  the  following  difficult  yet  realistic  situations  often  encountered  in  seismic  monitoring  for  the  purposes  of 
hazard  reduction  as  well  as  for  the  CTBT  verification:  1)  the  crustal  model  is  not  known,  and  hence  an  arbitrary  global 
model  (such  as  IASP91)  may  have  to  be  utilized;  2)  the  seismic  network  is  not  calibrated,  and  hence  the  phase  arrivals 
will  be  used  as  they  are  -  without  any  station-specific  or  path-specific  correction;  and  3)  the  azimuthal  coverage  of  the 
recording  stations  is  poor. 


Summary 

Two  non-standard,  independently  developed,  azi¬ 
muth-based  seismic  triangulation  procedures  are  im¬ 
plemented  and  evaluated  in  this  study. 

Through  azimuthal  triangulation,  Asian  seismic 
network  operators  in  China  and  lapan  have  been  locat¬ 
ing  earthquakes  for  several  decades  using  the  “Induced 
Perpendicular  Bisectors”  (IPB),  long  before  computers 
were  invented  and  introduced  to  seismologists.  The  un¬ 
derlying  simple  principle  can  be  best  illustrated  with  the 
following  extreme  case.  If  two  seismographs  happen  to 
record  identical  arrival  times  of  the  same  seismic  phase, 
then  to  the  extent  the  one-dimensional  (1-D)  structure 
applies  to  the  region,  the  hypocenter  should  lie  on  the 
perpendicular  bisector  of  the  line  segment  or  the  great 
circle  which  connects  these  two  specific  seismographs. 
Depending  on  the  epicentral  distances,  the  perpendicu¬ 
lar  bisector  itself  could  be  a  great  circle  along  the 
Earth's  surface  or  a  normal  section  cutting  through  the 
Earth.  If  two  or  more  such  perpendicular  bisectors  are 
available,  then  the  hypocenter  or  epicenter  can  be  de¬ 
termined  via  triangulation.  The  merit  lies  in  how  the 
perpendicular  bisector  is  derived  in  a  more  general  set¬ 
ting,  when  the  arrival  time  varies  from  station  to  station, 
which  is  typically  the  case  in  reality;  and  when  wave¬ 
form  data  are  not  available,  which  renders  techniques 
based  on  full  waveforms  (such  as  polarization  analysis 
and  frequency-wave  number  (FK)  technique)  not  appli¬ 
cable.  The  so-called  “Yin  Zhong  Xian”  (“3lct:,££”  in 
Chinese)  algorithm,  hereafter  the  YZX  method,  is  an 
Oriental  version  of  IPB-based  procedure.  It  computes 
one  IPB  (via  interpolation)  for  each  group  of  three 
seismographs. 

Also  relying  on  the  azimuthal  triangulation  for  seis¬ 
mic  location,  R.  Jih  [1]  proposes  a  procedure,  JO,  to  de¬ 
rive  the  back  azimuth  with  a  large -aperture  network 
where  all  seismographs  are  on  one  side  of  the  event. 
Any  standard  Geiger-type  of  least-squares  inversion 
routine  can  be  applied  to  determine  the  back  azimuth 
easily.  Two  or  more  such  skewed  networks  would  suf¬ 
fice  to  derive  two  back  azimuths  for  triangulation  pur¬ 
pose.  It  has  been  demonstrated  in  [1]  that  this  simple. 


hybrid  procedure  is  particularly  suitable  for  the  seismic 
location  problem  at  regional  distances  when  [2]  the 
crustal  model  is  not  known,  [3]  the  seismic  network  is 
not  calibrated,  and  [4]  the  azimuthal  coverage  of  record¬ 
ing  stations  is  poor.  In  this  paper  the  procedure  JO  is 
compared  against  the  YZX  method,  to  relocate  earth¬ 
quakes  and  explosions  of  known  or  well-constrained 
locations.  It  is  shown  that,  at  regional  distances,  the  JO 
algorithm  performs  better  than  does  the  YZX  method, 
because  the  back  azimuth  derived  with  JO  is  more  relia¬ 
ble  and  stable.  This  is  not  surprising,  as  the  YZX  tends 
to  be  more  susceptible  to  the  small-scale  lateral  hetero¬ 
geneity,  which  is  often  present  in  the  crust.  At  local  dis¬ 
tances,  however,  the  two  methods  seem  to  perform 
equally  well. 

1.  Background 

Epicenter  remains  as  one  of  a  few  important  source 
parameters  that  need  to  be  promptly  determined  with 
sufficient  accuracy  and  confidence  whenever  a  natural  - 
or  man-made  -  seismic  event  of  concern  occurs.  Over 
the  years  the  computer-based  Geiger-type  least-squares 
inversion  and  its  variations  have  become  the  standard 
means  in  solving  the  location  problem. 

If  the  ground-truth  location  of  a  nearby  seismic 
event  is  known,  the  arrival  readings  of  the  new  event  at 
each  seismograph  can  be  individually  “calibrated”  to 
match  the  predicted  travel-time  curves  of  some  pre¬ 
selected  velocity  model.  It  has  been  long  established 
that,  with  the  station  corrections,  locations  of  new  seis¬ 
mic  events  in  the  vicinity  around  the  ground-truth  event 
can  be  significantly  improved.  This  is  the  standard  ap¬ 
proach  of  “calibration”.  In  a  special  application  of  this 
technique,  travel-time  residuals  are  obtained  from  a 
least-square  location  routine  which  is  run  with  the  depth 
constrained  to  the  known  “true”  value.  When  these  re¬ 
siduals  are  used  as  travel-time  corrections  in  the  same 
program  run  depth-free,  nearby  events  can  be  located 
with  smaller  errors  in  depth.  An  elaboration  of  this  spe¬ 
cific  application  has  been  denoted  the  SRST  (Source- 
Region-Station-Time)  technique  by  [5]  (see  also  [6]). 
Alternatively,  a  postulated  ID  crustal  velocity  model 
can  be  adjusted  so  that  the  estimated  hypocenter  would 
perfectly  match  the  ground-truth  event.  Thus,  the  fine- 
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tuned  velocity  model  can  be  used  to  determine  the  loca¬ 
tion  of  nearby  events  without  the  need  of  establishing 
residual  corrections  for  individual  stations.  Both  the 
SRST  and  the  model-tuning  procedures  require  at  least 
one  ground-truth  event.  The  North  Korean-declared 
nuclear  tests  at  Punggye  serve  as  a  classic  example 
where  multiple-event  location  techniques  would  work 
very  well,  and  the  impressive  seismic  results  seem  to  be 
corroborated  by  satellite  interferomeUy  [2], 

There  are  cases  where  no  nearby  event  of  known  lo¬ 
cation  is  available.  Furthermore,  there  are  cases  where 
not  only  the  crustal  model  is  not  well  known,  the  seis¬ 
mic  stations  could  be  sparse  or  even  spread  in  a  very 
skew  geometry.  R.  Jih  [1]  presents  a  procedure,  JO, 
which  could  be  particularly  suitable  for  these  disadvan¬ 
tageous  situations.  In  this  study,  JO  and  another  azi¬ 
muth-based  non-standard  location  procedure  are 
reviewed  and  compared. 

2.  Adaptive  Location  Methods  JO  and  YZX 

The  most  commonly  used  triangulation  procedure 
for  location  is  the  “radial  triangulation”  which  is  based 
on  intersection  of  several  arcs  drawn  at  respective  epi- 
cenUal  distances.  The  epicentral  distances  can  be  indi¬ 
vidually  estimated  at  each  seismograph  if,  for  instance, 

[2]  the  differential  arrival  (such  as  S-P)  is  available  and 

[3]  the  P-  and  S-wave  velocity  structure  around  the 
seismograph  is  known.  In  the  case  only  one  common 
phase,  say,  the  initial  P,  is  available,  then  each  pair  of 
seismographs  can  jointly  determine  a  hyperbola  which 
passes  through  the  epicenter  -  provided  that,  again,  the 
P-wave  velocity  is  known.  Lacking  the  knowledge  of 
the  crustal  structure,  the  location  derived  from  radial 
triangulation  could  be  severely  biased  when  the  number 
of  seismograph  is  limited.  The  Geiger-type  least  square 
inversion  using  arrival  times  would  inherently  suffer 
from  the  same  drawback  and  limitations. 

Through  “azimuthal  triangulation”,  Asian  seismic 
network  operators  have  been  locating  earthquakes  for 
several  decades  using  the  “Induced  Perpendicular  Bisec¬ 
tors”  (IPB),  even  before  computers  were  introduced  to 
seismologists.  The  underlying  simple  principle  can  be 
best  illustrated  with  the  following  extreme  case.  If  two 
seismographs  happen  to  record  identical  arrival  times  of 
the  same  seismic  phase,  then  to  the  extent  the  ID  sttuc- 
ture  applies  to  the  region,  the  hypocenter  should  lie  on 
the  perpendicular  bisector  of  the  line  segment  or  the 
great  circle  which  connects  these  two  specific  seismo¬ 
graphs.  Depending  on  the  epicentral  distances,  the  per¬ 
pendicular  bisector  itself  could  be  a  great  circle  along 
the  Earth's  surface  or  a  normal  section  cutting  through 
the  Earth.  If  two  or  more  such  perpendicular  bisectors 
are  available,  then  the  hypocenter  or  epicenter  can  be 
determined  via  triangulation.  The  merit  lies  in  how  the 
perpendicular  bisector  is  derived  in  a  more  general  set¬ 
ting  when  the  arrival  time  varies  from  station  to  station, 
which  is  typically  the  case  in  reality;  and  when  wave¬ 
form  data  are  not  available,  which  renders  techniques 
based  on  full  waveforms  (such  as  polarization  analysis 


and  frequency-wavenumber  (FK)  technique)  not  appli¬ 
cable. 

The  so-called  “Yin  Zhong Xian’’  algorithm,  hereafter 
the  YZX  method,  is  an  oriental  version  of  IPB-based 
procedure.  Given  three  arrival  picks,  tl  <  t2  <  t3,  at 
three  seismographs  SI,  S2,  and  S3,  respectively,  the 
YZX  method  finds  a  point  S4  which  lies  on  the  great 
circle  connecting  SI  and  S3  such  that  the  arrival  time  at 
S4  would  be  approximately  t2.  S2  and  S4  then  deter¬ 
mine  one  IPB.  The  search  of  S4  (through  interpolation) 
requires  a  priori  knowledge  about  the  velocity  structure. 
In  reality,  simple  interpolation  procedures  are  often 
used,  which  inevitably  introduce  additional  uncertainty 
into  the  derived  back  azimuth,  and  hence  the  location 
would  be  affected  as  well. 

Also  relying  on  azimuthal  triangulation  for  seismic 
location,  R.Jih  [1]  proposes  a  procedure  JO  to  derive  the 
back  azimuth  with  a  large  aperture  network  where  all 
seismographs  are  on  one  side  of  the  event.  The  basic 
idea  of  JO  method  is  to  decompose  the  recording  net¬ 
work  into  (at  least)  two  sub-networks  so  that  each  sub¬ 
network  is  comprised  of  stations  spread  over  only  one 
side  of  the  postulated  epicenter.  That  is,  each  sub¬ 
network  should  have  a  very  large  azimuthal  gap  and 
preferably  a  large  aperture  as  well.  Any  standard  Gei- 
ger-type  of  least-squares  inversion  routine  can  be  ap¬ 
plied  to  determine  the  back  azimuth  easily.  Two  or  more 
such  skew  networks  would  suffice  to  derive  two  back 
azimuths  for  triangulation. 

In  this  paper  the  procedure  JO  is  compared  against 
the  YZX  method,  to  relocate  earthquakes  and  explo¬ 
sions  of  known  or  well-constrained  locations.  It  is 
demonstrated  that,  at  regional  distances,  the  JO  algo¬ 
rithm  performs  better  than  does  the  YZX  method, 
because  the  back  azimuth  derived  with  JO  is  more 
reliable  and  stable.  This  is  not  surprising,  as  the  YZX 
tends  to  be  more  susceptible  to  the  small-scale  lateral 
heterogeneity  which  is  often  present  in  the  crust.  At 
local  distances,  however,  the  two  methods  seem  to  per¬ 
form  equally  well  (not  shown). 

3.  Example  1:  The  Salmon  Explosion 

of  October  22, 1964 

Salmon,  an  underground  nuclear  explosion  in  Mis¬ 
sissippi,  was  seismically  recorded  throughout  North 
America  and  at  some  teleseismic  stations.  The  5-kt  ex¬ 
plosion  was  detonated  at  16h  00m  00.0s  UT  on  October 
22,  1964,  at  a  depth  of  828.1  meters.  The  epicenter  is 
31°08’31.57”N,  89°34’11.8”W,  in  the  Taut  salt  dome  in 
Hattiesburg,  Southern  Mississippi.  Extensive  studies 
have  been  conducted  for  this  event  to  validate  crustal 
profiles  toward  different  directions  [7  -  9,  J.  Jordan  [et 
al.]  in  [7]  published  phase  picks  recorded  at  143  stations 
for  Salmon  explosion.  Thirty  six  arrivals  from  those 
stations  within  10°  are  used  in  this  example. 

Figure  1  shows  21  “induced  bisectors”,  each  derived 
with  the  YZX  method  from  three  sites  of  the  36-station 
regional  network.  Figure  2  is  the  same  as  Figure  1,  ex¬ 
cept  that  25  back  azimuths  are  determined  using  the 
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method  JO  in  conjunction  with  a  global  continental  av¬ 
erage  model,  IASP91  [10].  The  IASP91  model  is  arbi¬ 
trarily  chosen  to  mimic  the  situation  where  the  true 
crustal  structure  of  a  region  of  interest  may  not  be  avail¬ 
able,  and  seismic  analysts  would  have  to  use  a  possibly 
erroneous  crustal  model  to  start  with  (Figure  2).  The 


G.T.  location=31.142N -89.570W 0.828 
21  “induced  bisectors”  used,  each  derived  from  a  3-station  subnetwork 


method  JO  gives  tightly  converged  rays,  with  the  majori¬ 
ty  intersecting  in  the  vicinity  of  the  ground-truth  loca¬ 
tion.  The  dashed  circle  in  Figures  1,  2,  3,  4,  and  6  has  a 
radius  of  18.73  km,  which  encircles  an  area  of  1,000 
km2. 


-91.750  -67.390 


G.T.  location=31  142N  -89.570W  0.828km 
IASP91  model  used  for  all  subnetworks 


Figure  1.  Adaptive  Location  of  Salmon  w.  Method  YZX 

4.  Example  2:  The  Taiwan  Strait  Earthquake 
of  September  16, 1994 

The  September  16,  1994,  Taiwan  Strait  earthquake 
caused  significant  damages  to  southeast  China.  Though 
at  equal  distance  to  Taiwan,  its  damage  to  Taiwan  was 
relatively  minor.  Paths  to  the  mainland  China  apparently 
had  a  weaker  attenuation  in  ground  motion.  Fifteen 
pairs  of  Pn  and  Sn  arrival  times  have  been  published  by 
the  Fujian  Provincial  Seismological  Bureau  [FSB],  Chi¬ 
na  [11,  12],  W-Y.  Yeh  [11]  selected  five  sets  of  IPBs 
and  suggested  that  YZX  would  lead  to  a  solution  con¬ 
sistent  with  the  conventional  triangulation  results  using 
solely  S-P  times  (c.f.  Figures  2  and  3  of  [11]).  In  an 
attempt  to  duplicate  Yeh's  result,  it  is  found  that  the 
YZX  method  degrades  somewhat  as  additional  IPBs  are 
drawn  from  different  combinations  of  FSB  seismo¬ 
graphs  (Figure  3).  It  should  be  noted  that,  however,  the 
station  coordinates  of  FSB  seismographs  used  in  this  re¬ 
location  exercise  are  not  exact,  as  they  were  digitized 
from  geographic  maps,  which  certainly  have  added  ex¬ 
tra  random  errors  to  the  location,  with  an  effect  equiva¬ 
lent  to  those  of  larger  phase  picking  errors  or  stronger 
lateral  heterogeneities  in  the  crust. 

If  the  same  1 -dimensional  model  is  assumed  for 


Figure  2.  Adaptive  Location  of  Salmon  w.  Method  JO 

paths  to  Taiwan  and  mainland  China,  then  combining 
FSB  phase  picks  with  those  measured  at  Taiwan's  Cen¬ 
tral  Weather  Bureau  [CWB]  leads  to  a  location  with 
large  offset  (See  the  yellow  diamond  in  Figure  4).  The 
1 -second  difference  between  FSB  and  CWB  clocks,  as 
speculated  by  [12],  cannot  count  for  this  dramatic  shift 
of  location.  It  appears  that  paths  to  Fujian  region  may 
have  a  velocity  faster  than  that  of  the  IASP91  model, 
whereas  paths  to  Taiwan  may  be  much  slower.  This  is 
evidenced  by  Figure  4  (bottom)  in  which  FSB  and  CWB 
underestimate  and  overestimate,  respectively,  the  epi- 
central  distances  when  the  IASP91  model  is  used.  Nev¬ 
ertheless,  the  individual  back  azimuths  determined  by 
FSB  and  CWB  networks  are  not  severely  biased  in  di¬ 
rection.  Thus  applying  the  adaptive  method  JO  would 
lead  to  a  triangulation  result  extremely  close  to  the  ISC 
and  USGS/NEIC  solution,  which  are  based  on  hundreds 
of  stations  spread  over  all  directions.  The  adaptive 
method  is  less  sensitive  to  the  choice  of  crustal  model. 
Changing  the  crustal  model  from  IASP91  to  the  south¬ 
east  China  model,  which  is  in  routine  use  at  FSB  [13] 
(see  also  [14]),  does  not  seem  to  affect  the  triangulation 
result  (Table).  This  is  yet  another  indication  that  the 
adaptive  method  JO  is  robust. 
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Table.  Comparison  of  Epicenters  of  September  16,  1994  Event 


Network  /  Bulletin 

<P°(N) 

r  (E) 

Remark 

Fujian  Network 

22.60 

118.68 

[1 1],  Southeast  China  model 

Fujian  +  Guangdong 

22.70 

118.75 

[11],  Southeast  China  model 

CSB  Preliminary 

22.60 

118.73 

(China  Seismological  Bureau) 

CSB  RRSN 

23.00 

118.50 

(CSB  Rapid  Reporting  Seismograph  Network) 

Taiwan  CWB 

22.43 

118.47 

(Taiwan  Central  Weather  Bureau) 

IESAS  +  NCU 

22.37 

118.63 

[12] 

USGS/NEIC 

22.53 

118.71 

ISC 

22.52 

118.75 

Adaptive  Method  JO 

22.48 

118.73 

This  study:  Fujian+Taiwan,  SE  China  model 

Adaptive  Method  JO 

22.48 

118.72 

This  study:  Fujian+Taiwan,  IASP91  model 

(Note:  China  Seismological  Bureau,  or  CSB,  was  the  name  of  the  Chinese  earthquake  monitoring  agency  at  the  time  the  earthquake  occurred, 
which  has  since  been  renamed  China  Earthquake  Administration,  or  CEA.) 


30.000 


21.000 


118.000  119.250 


NEIC=0620  18.7  22.53N  118.71E  13km,  ISC=22.52N  118.75E  13  km; 
CWB=0620  15.62  22.43N  118.47E  19km;  FSB=22.60N  118.68E  10  km; 
10  “induced  bisectors”  used,  each  derived  from  a  3-station  subnetwork 


106.000  123.000 


118.000  110.260 


NEIC=0620  18.7  22.53N  118.71E  13km;  ISC  (REF)=22.52N  118.75E  13  km; 
CWB  (Taiwan)=0620  15.62  22.43N  118.47E  19,1  km;  FSB=0620  22.10.22.60N 
118.68E  10,0  km;  Adaptive  JO  w.  IASP91: 22.48N  118.72E 


Figure  3.  Adaptive  Location  of  940916  w.  Method  YZX 


Figure  4.  Adaptive  Location  of  940916  w.  Method  JO 


5.  Example  3:  The  Kara  Sea  Event 
of  August  16, 1997 

Figure  5  illustrates  how  the  method  JO  can  be  ap¬ 
plied  to  derive  the  location  of  Kara  Sea  event  of  August 
16,  1997.  A  sub-network  “A”  comprised  of  6  stations 
(KBS,  SPITS,  KEV,  SDF,  KAF,  and  FINESS),  11  phase 
picks,  are  fed  into  the  LocSAT  program  (Figure  5,  top) 
as  described  in  [15,  16].  The  semi-major  axes  of  the 
seven  error  ellipses  are  nearly  the  same,  all  pointing  to 
the  east  -  roughly  100°  from  the  north  (Figure  6).  A 
second  sub-network  "  B"  is  formed  with  7  stations  all 
spread  to  the  south  of  Novaya  Zemlya,  which  includes 


NRI,  ARU,  PKK,  JOF,  VAF,  SUF,  and  KJN.  There  are 
10  phase  picks  reported  by  these  stations.  Seven  differ¬ 
ent  continental  crustal  models  produce  nearly  identical 
“optimal”  back  azimuths,  pointing  to  the  north  (Figure 
5,  bottom).  Using  the  IASP91  model,  the  two  rays  inter¬ 
sect  at  a  location  1 1  km  to  the  west  of  pIDC's  REB  loca¬ 
tion  published  in  [19]  (Figure  6).  If  the  two  networks 
were  combined  and  the  conventional  location  algorithm 
applied,  the  resulting  location  is  24  km  to  the  south  of 
the  REB  location  (unfilled  star  in  Figure  6).  In-depth 
reviews  addressing  both  the  technical  and  political  as¬ 
pects  of  this  event  can  be  found  in  [18,  19]. 
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Network  A=FINES  KAF  KBS  KEV  SDF  SPITS 
Network  B=ARU  JOF  KJN  NRI  PKK  SUF  VAF 


Figure  5.  Grouping  of  Seismographs  for  970816 


]REB  location  =  72.6484  57.3517  [depth  set  to  0]  (Cross) 

IASP91, 11  picks  ==>  OT=2:1 1:5.281,  72.750  55.358  0.000,  Offset=67.2 
(Network  A=FINES  KAF  KBS  KEV  SDF  SPITS) 

IASP91, 10  picks  ==>  OT=2:1 1:4.357, 71.839  57.472  55.559,  Offset=90.4 
(Network  B=ARU  JOF  KJN  NRI  PKK  SUF  VAF) 

IASP91  +  JO  method  =>  72.657  57.033,  Offset=10.7  (Blue  Star) 

IASP91  +  all  picks  combined  =>  72.445  57.104,  Offset=24.2  (Blank  Star) 

Figure  6.  Locations  of  970816,  Conventional  vs.  Adaptive  [JO] 


6.  Discussion  and  Conclusions 

Two  azimuth-based  location  techniques,  JO  and 
YZX,  have  been  implemented  and  tested  to  relocate 
earthquakes  and  explosions  of  known  or  well- 
constrained  locations.  It  is  shown  that,  at  regional  dis¬ 
tances,  the  JO  algorithm  performs  better  than  does  the 
YZX  method,  because  the  back  azimuth  derived  with  JO 
is  more  reliable  and  stable.  This  is  not  surprising,  as  the 
YZX  tends  to  be  more  susceptible  to  the  small-scale 
lateral  heterogeneities  which  are  often  present  in  the 
crust.  At  local  distances,  however,  the  two  methods 
seem  to  perform  equally  well. 

An  interesting  question  raised  by  Prof.  Thorne  Lay 
at  the  1st  Event  Screening  Workshop  (November  1997, 
Beijing)  was  that,  if  both  the  staged  location  procedure 
(like  JO)  and  the  conventional  Geiger  inversion  are  us¬ 
ing  virtually  the  same  amount  of  information,  would  the 
staged  approach  offer  advantages  over  the  conventional 
procedure?  In  the  case  of  simple  crustal  structure,  in¬ 
deed  both  methods  could  perform  equally.  Figure  2  ex¬ 
hibits  a  tight  clustering  of  back  azimuths  derived  from 
the  adaptive  method  JO,  resulting  a  well-constrained 
epicenter.  This  can  certainly  be  achieved  with  all  the 
phase  picks  lumped  in  one  single  Geiger  inversion.  This 
is  not  the  case  in  general,  however.  As  shown  in  Figure 
4,  lumping  phase  picks  from  Fujian  and  Taiwan  seismic 
networks  into  one  least-squares  inversion  would  lead  to 
a  disastrous  location,  unless  the  difference  in  lateral 
difference  in  crustal  velocities  is  accounted  for.  This  is 
the  intrinsic  drawback  of  the  conventional  location  pro¬ 
cedure.  In  the  case  of  staged  procedure  such  as  JO,  each 
back  azimuth  is  determined  with  a  skew  sub-network. 
The  skew  geometry  of  the  sub-network  would  inevita¬ 
bly  lead  to  large  uncertainty  in  location  and  large  error 
in  epicenter.  However,  the  method  JO  does  not  use  the 
epicenter  information  in  the  second  stage  of  the  process. 
The  only  piece  of  information  -  which  is  the  best  con¬ 
strained  piece  -  utilized  by  JO  is  the  semi-major  axis,  or 
equivalently,  the  back  azimuth.  As  a  result,  the  JO  meth¬ 
od  is  practically  using  several  back  azimuths  -  each  one 
being  better  constrained  in  the  process  -  for  triangula¬ 
tion.  This  explains  why  JO  is  not  equivalent  to  the  con¬ 
ventional  Geiger  inversion. 

A  side  remark  on  seismic  data  should  be  made. 
Seismic  data  have  multiple  uses  -  quick  and  accurate 
location  of  damaging  earthquake  so  as  to  permit  rapid 
emergency  response,  studies  of  the  interior  of  the  earth 
and  the  physics  of  the  earthquake  source,  and  nuclear 
test  verification  [18].  The  field  of  seismology,  which 
has  been  involved  in  earthquake  studies  for  nearly  100 
years,  has  a  long  tradition  of  international  data  ex¬ 
change,  much  like  that  for  weather  information.  It 
would  be  a  tragedy  if  a  major  earthquake  disaster  strikes 
and  response  to  the  disaster  were  delayed  because  seis¬ 
mic  data  were  not  available.  Example  2  clearly  illus¬ 
trates  the  advantages  and  importance  of  sharing  the 
phase  readings  between  China  and  Taiwan.  China 
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Earthquake  Administration  [CEA,  formerly  China 
Seismological  Bureau  or  CSB]  has  been  promptly  shar¬ 
ing  phase  readings  of  large  earthquakes  with  Russia, 

Pakistan,  Algeria,  Laos,  Myanmar,  Samoa,  Indonesia, 
and  Cuba  [2.4],  It  is  not  surprising  that  the  performance 
of  CEA's  RRSN  [Rapid  Reporting  Seismograph  Net¬ 
work]  in  earthquake  location  has  greatly  benefited  from 
the  data  exchange.  News  has  it  that  CEA  and  the  Cen¬ 
tral  Weather  Bureau  [CWB]  of  Taiwan  have  signed  a 
seismic  cooperation  agreement  in  February  2014,  and 
have  started  routine  seismic  data  exchange  between  the 
two  sides  of  the  Taiwan  Strait  since  June  2015. 

A  lesson  learned  from  Example  2  (Figure  4)  is  that 
phase  picks  from  different  tectonic  regions  must  be 
combined  in  a  judicious  manner,  if  the  conventional 
Geiger  inversion  is  to  be  applied.  One  way  to  do  so  is  to 
take  into  account  the  difference  in  crustal  structures  and 
apply  the  path-dependent  travel-time  tables  -  or  equiva¬ 
lently,  path-dependent  corrections  -  prior  to  lumping  the 
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phase  picks  together  for  inversion.  If  adaptive  methods 
(such  as  the  method  JO)  are  utilized,  however,  then  the 
improvement  in  location  can  be  immediately  achieved 
without  the  need  to  wait  for  the  time-consuming  "cali¬ 
bration  research"  is  conducted.  This  can  be  regarded  as 
an  ad  hoc  measure. 
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A3HMYTTM  ECEIITEy  HEri  JIJW  O^HFAJIAPflbl  CTAHflAPTTbl 
EMECE  ^CEPriJIIKTEy^m  EKI  0flICTEPIH  CAJIblCTblPMA  3EPTTEY 

jJ'/KH  P. 

AKJU  MeMjieKemmiK  denupmuMenminin  ljapyjiany  niypanbi  tcejiicwdepdi  oapu.iay, 
eepu (/> m<an uiijiay  .wane  caipnay  6wpocbi,  Bauiummon,  AKJU 

VtaKana/ia,  atitNiyTTbi  ccenTey  iliac  nern/ieareir,  etci  CTair/iapTTbi  eMec,  6ip  6ipme  oaiiaaiibiccbi'i  33ipaeHreH 
ceiJcMHKajiBiK;  TpnaHrynaunji  npoueziypacti  KapacTBiptuia/iBi.  3epTTeynepz(iH  Mascara,  KayinTi  TeMeHxieTy,  coH^aH-ai? 
JIOKTIII  caKTayuit  oaKbiaay  vtaKcaTbin;iaFbi  cciicMitKaabrK  vtonuTop itnricinac  Keneci  naKTbi,  mti  KCi.acccTin 
yKarpafijiappa  OKHFajiap,nbi  ceitcMHicajibn?  >KepriniKTey  npo6jieMajiapi>iH  TuiM^i  meiny  yuim  anropiiTM^i  33ipneyi 
bonraH:  1)  *ep  KbipTbicbiHbm  Modern  6eHM3:iiM,  conntncraH,  epKiHine  ritobajib  Mo^eni  KOn/taHtuia^bi  (Mbicajibi, 
IASP91);  2)  CTaiouMJiapxtbm  cencMHicajibiK;  >Keiiici  KanH6pneH6ereH,  coH/tbnyraH,  tj)a3anap  Tycy  yaKbrrrapbi  6acTam?i>i 
TypiHzte  KOn/iaHi>raa/ii>i  -  CTaHUHtuibiK  Ty3eTy:iepci3  >K3He  5Kyripy  >K0itbiHa  Ty3eTy:iepci3;  3)  TipKeimH  CTaHinuiaapMeH 
o^HFaitap^bi  a3HMyTaitztiK  Kaivnybi  Haniap. 

CPABHHTEJIbHOE  HCCJIE^OBAHHE  flBYX  METOflOB  HECTAH^AI’THOH 
JIOKAJIH3AIIHH  COBbITHH  HA  OCHOBE  PACHETA  A3HMYTA 

JJ/KII  P. 

Ewpo  KomnpojiH,  eepwpiiKanuu  u  codnmdenua  cozjiaiuenuu  o  ttoopy.nceniiiix 
rocydapctneeHuozo  denapmaMenma  CUIA,  Buuuinzmou,  CU1A 

B  CTaTbc  p a c c m a t p it b a i o t c h  /(Be  irccTaii/iapTiibic,  paipaooTaniibic  He3aBHCHMO  ztpyr  ot  apyra,  npoiteaypbi 
cciicMit'iccKoii  Tpiiaiiryaaiuiii.  ocnoBainibic  Ha  pac'iCTC  a3HMyTa.  Ijcabio  itccac/io  Bailin')  Obiaa  pa3pa6oTKa  anropHTMa 
pirn  3tj)(|)eKTHBHoro  pemeHHJi  npo6aeMbi  ceiicMHHecKOH  jioKajiH3auHH  co6mthh  b  cite/iyioinHx  peajiBHbix,  uacro 
BCTpC'iaiOIHHXCSt  CHTyaUHStX  npn  CCitCMH'ICCKOM  MOHHTOpHHre  RJIfL  CHH)KeHHJI  OnaCHOCTH,  a  T3K)Ke  pjm  KOHTpOJIH  3a 
co6aio/iciiiic\t  /I  B  3  JIM:  1)  Koponaa  \to/icab  He  H3BecTHa,  cacyoBaTcabiio.  npiiMcnstcTcst  npoittBoabiiaa  raooaabitaa 
Mozieat  (HanpiiMep,  IASP91);  2)  cciicMii'iccKaa  ccTb  cTaniutii  He  OTKaaitoponana,  cac.aoBaTcabno,  BpeMeHa 
BCTynaenitst  (fiat  itcnoabtyiOTca  b  hx  itcxoanovt  Bit.ac  -  6e3  CTanmioimbix  nonpaBOK  h  6e3  nonpaBOK  Ha  nyTb  npo6era; 
3)  aiitMyTaabiibiii  oxBaT  coObiTitst.  pcntcTpiipyioiuiiMit  CTaimiistMit  naoxoii. 


89 


nEPI/IOflMHECKMM  HAyMHO-TEXHMHECKMM  XCyPHAJl 
HAUMOHAflbHOrO  flflEPHOI~0  UEHTPA 
PECnYEJlMKM  KA3AXCTAH 


ISSN  1729-7516 


